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Abstract. The effect of aerosol on clouds poses one of the 
largest uncertainties in estimating the anthropogenic contri- 
bution to climate change. Small human-induced perturba- 
tions to cloud characteristics via aerosol pathways can create 
a change in the top-of-atmosphere radiative forcing of hun- 
dreds ofWm" 2 . Here we focus on links between aerosol and 
deep convective clouds of the Atlantic and Pacific Intertrop- 
ical Convergence Zones, noting that the aerosol environment 
in each region is entirely different. The tops of these ver- 
tically developed clouds consisting of mostly ice can reach 
high levels of the atmosphere, overshooting the lower strato- 
sphere and reaching altitudes greater than 16 km. We show 
a link between aerosol, clouds and the free atmosphere wind 
profile that can change the magnitude and sign of the overall 
climate radiative forcing. 

We find that increased aerosol loading is associated with 
taller cloud towers and anvils. The taller clouds reach levels 
of enhanced wind speeds that act to spread and thin the anvil 
clouds, increasing areal coverage and decreasing cloud opti- 
cal depth. The radiative effect of this transition is to create a 
positive radiative forcing (wanning) at top-of-atmosphere. 

Furthermore we introduce the cloud optical depth (r), 
cloud height (Z) forcing space and show that underestima- 
tion of radiative forcing is likely to occur in cases of non 
homogenous clouds. Specifically, the mean radiative forcing 
of towers and anvils in the same scene can be several times 
greater than simply calculating the forcing from the mean 
cloud optical depth in the scene. 
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Limitations of the method are discussed, alternative 
sources of aerosol loading are tested and meteorological vari- 
ance is restricted, but the trend of taller clouds, increased and 
thimier anvils associated with increased aerosol loading re- 
mains robust through all the different tests and perturbations. 


1 Introduction 

Increase in aerosol concentration may invigorate convective 
clouds through the particles’ role as cloud condensation nu- 
clei (CCN) and the related feedbacks ignited by changes of 
the clouds 4 initial droplet size distribution (Andreae et al., 
2004; Koren et al., 2005). The proposed chain of events links 
initial changes in the size distribution of the cloud’s droplets 
to changes m the net updraft velocities, droplet growth rates, 
cloud vertical development and precipitation in the warm, 
mixed and cold phases. 

Previous studies show clear correlations between cloud 
coverage, lifetime (Albrecht, 1989; Jacobson et al., 2007; 
Kaufman and Koren, 2006) and optical properties (Twomey, 
1977; Rosenfeld, 2000). It was shown from satellite data 
analysis and in situ measurements, that for convective clouds, 

increase m aerosol loading is associated with taller invig- 
orated clouds, larger cloud fraction and more extensive ice 
portions. These associations are found over the tropical At- 
lantic (Koren et al., 2005; Jenkins and Pratt, 2008; Jenkins et 
al., 2008), Europe (Devasthale et al., 2005), North and South 
America (Andreae et al., 2004; Koren et al., 2008a; Lind- 
sey and Fromm, 2008), and appear for all types of aerosol 
particles: biomass burning smoke, urban/industrial aerosol 
and desert dust. Numerical modeling studies show aerosols 
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penetrating cloud base of maritime clouds dramatically in- 
creasing the amount of supercooled water, as well as the 
ice contents and vertical velocities (Khain et al., 2008a, b). 
Aerosol particles can also affect convective cloud properties 
when ice nucleating particles are entrained directly into the 
ice portion of the cloud (Fridlind et al., 2004). 

Deep convective clouds consist of the cloud tower (or tow- 
ers) where most of the convection occurs and the cloud anvil 
(or anvils), which is a thinner layer of mostly ice particles ex- 
tending away from the cloud towers. Anvils form when part 
of the cloud spreads out below a stable atmospheric (inver- 
sion) layer, such as the tropopause. The anvil horizontal area 
can be much larger than the tower area, as the inversion layer 
above it acts as a ceiling, blocking the anvil cloud elements 
from expanding upwards and forcing horizontal expansion. 
The anvil elements are further spread by the horizontal winds 
below the inversion level (Fig. la). 

Over the tropical Atlantic, convective cloud liquid water 
optical depth and ice optical depth respond differently to 
increase in aerosol loading (measured by the aerosol opti- 
cal depth, AOD) as observed from satellites (Koren et al., 
2005). As expected, due to changes in the droplet size- 
distribution and warm rain suppression, liquid cloud opti- 
cal depth increases with AOD, however, ice cloud optical 
depth decreases. Here by analyzing horizontal wind pro- 
files from radiosonde data from several stations in the tropics 
and subtropics (NOAA/ESRL Radiosonde Database), cloud 
and aerosol level-3 MODIS version 5, daily data (MODIS - 
Moderate Resolution Imaging Spectroradiometer, Platnick et 
al., 2003; Remer et al., 2008), we show that the clear reduc- 
tion in the average ice optical depth with increase in aerosol 
loading is due to a strong increase in the ice anvil area. We 
then estimate the anvil enhancement consequences to the cli- 
mate radiative balance using a tropical profile (Ricchiazzi et 
al., 1998). The study area focuses on the tropical Atlantic 
(0° N to 14° N; 1 8° W to 45° W) and tropical Pacific (0° N to 
14° N; 180° W to 150° W) of June- July-August 2007 using 
the MODIS-Aqua data (01:30 p.m. local-tune overpass). 

2 Analysis and results 
2.1 MODIS data 

We introduce the Tower to Anvil area ratio (TAR), and show 
that while clouds are invigorated by increasing aerosol load- 
ing and reach higher levels in the atmosphere, this ratio 
decreases due to the associated increase and spread of the 
anvil area. The anvil area increases both due to the en- 
hanced convection and stronger horizontal winds in the up- 
per troposphere. Figure lb shows a whole month average 
(August) daytime horizontal wind speed of the free atmo- 
sphere (above the boundary layer) as a function of pressure 
level for 3 stations along the tropical Atlantic: Natal, Brazil 
(5.9° S 35.2° W): Illia Fernando de Noronha, Brazil (3.8° S 


32.4° W), Ascension Island (7.9° S 14.4° W). A consistent 
sharp increase with height occurs from 300 to 150 hPa (~9 to 
15 km) and then decreases towards the tropopause (100 hPa). 
This behavior is consistent for the wide variety of stations 
analyzed (day and night, land and ocean, inside and outside 
the tropics, for all seasons). The high level jet contributes to 
spreadmg and diluting the anvils, expanding the total cloud 
cover and decreasing the TAR. 

The link between cloud fraction and cloud height (mea- 
sured as cloud top pressure) over the Atlantic study area is 
shown in Fig. lc. Indeed it is shown that the average cloud 
fraction (each point is an average of 60 one-degree MODIS 
pixels after sorting all pixels by their cloud top pressure val- 
ues) increases as a function of the cloud top height. The 
change in slope at ~850hPa marks the transition from the 
lower (marine boundary layer) to the free - upper tropo- 
sphere where the clouds have larger average fraction. At 
300 hPa almost complete overcast conditions are reached. 

How does aerosol loading affect the dependence of cloud 
fraction and height? After filtermg out shallow clouds (cloud 
tops below 800 hPa), the cloud fraction data over the Atlantic 
was sorted according to the aerosol optical depth (AOD) and 
divided into 3 equally-sampled groups with mean AOD of 
each group equal to 0.13, 0.24, 0.41, respectively. Then 
the same analysis described for Fig. lc was done for each 
subgroup (see for more details: Koren et al., 2005, 2008a). 
While the cloud fraction of all 3 groups increases as a func- 
tion of the pressure level, a clear shift of the polluted clouds 
towards larger fraction is shown for each pressure level. 
Moreover the invigoration effect is also shown as the pol- 
luted clouds are taller (Fig. Id). 

MODIS retrieves detailed cloud and aerosol properties 
with resolution of 1 to 10 km. The data are averaged into 
a daily 1° resolution grid allowing for simultaneous obser- 
vations of aerosols in cloud-free regions and clouds in the 
cloudy regions of the grid box when not completely overcast. 
To save the information of the higher resolution retrievals 
(1 km 2 for clouds and 10 km 2 for aerosol) histograms of the 
original cloud and aerosol properties are provided (King et 
al., 2003). The histograms of the ice cloud optical depth (r) 
were analyzed as a function of the AOD over the study ar- 
eas. After inspection of many high resolution anvil optical 
depth data over the region and based on the cloud optical 
depth distribution (Fig. 2), we divided the histograms into 
two regimes: anvil (r < 10) and tower (r > 10). The Tower 
to Anvil area ratio (TAR) is the averaged number of pixels in 
the tower regime ( N t ) divided by the number of pixels in the 
anvil regime (N a ): 


The average cloud optical depth z c is weighted by the hori- 
zontal coverage area of the anvil optical depth ( r a ) and the 
tower optical depth (r r ): 
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Fig. 1. (a) Photograph of a deep convective cloud with large anvil over Africa - 5 February 2008, photographed from the space shuttle (Image 
courtesy of the Image Science & Analysis Laboratory, NASA Johnson Space Center), (b) Monthly average daytime wind profiles measured 
with radiosondes at 3 stations along the tropical Atlantic: Natal, Brazil (red, 5.9° S 35.2° W); Ilha Fernando de Noronha, Brazil (black, 
3.8° S 32.4° W), Ascension Island (blue, 7.9° S 14.4° W). (c) Analysis of cloud fr action as a function of the cloud top pressure (height) over 
the ITCZ (tropical Atlantic), (d) Focusing on the free atmosphere (>800hPa), cloud fraction as a function of the cloud top pressure for 3 
averaged aerosol groups, AOD=0. 13 (blue), 0.24 (green) and 0.41 (red), for the tropical Atlantic. Note that the pressure level ranges (height) 
are different for b, c and d. 


N t T t +N a T a TAR- z, + x a 
N t + N a TAR+ 1 v ’ 

As the TAR becomes smaller, x a contributes more to the to- 
tal cloud optical depth r c . Figure 3a shows that the TAR 
decreases significantly (from ~0.6 to ~0.35) for an average 
AOD variation from clean M).l to polluted 0.3. Further- 
more, Fig. 3b shows that the average ice cloud fraction in- 
creases significantly as a function of the AOD and that most 
of the increase is in the anvil area. Figure 3c shows that as the 
aerosol loading increases, the anvil height increases and the 
anvil optical depth decreases significantly. The same analy- 
ses are shown for the Pacific dataset in Fig. 3d to f, but for a 
smaller AOD range (0.05 > AOD >0.2) more appropriate for 
the less aerosol-laden Pacific. Over the Atlantic we repeated 
the analysis using 2002 MODIS Terra (10:30 a.m. local-time 
overpass) and 2005 Aqua data. The trends were similar in all 
years both for Terra and Aqua. 

The proposed chain of events suggested here is the fol- 
lowing: aerosols serving as CCNs invigorate convection and 
increase the cloud vertical extent into the upper troposphere 
(Fig. Id, 3c and 3f). The taller clouds now reach higher in 
the atmosphere and spread horizontally, aided by the stronger 
winds aloft (Fig. lb), and result in higher (mostly anvil) 
cloud fraction (Figs. Id, 3b and e). As the wind shear 



Fig. 2. The histogram of the ice cloud optical depth (r) values over 
the Atlantic study area. The red dots mark the MODIS histogram 
bins. A r range between 6 to 18 was tested as a threshold between 
anvils to towers and all threshold showed the same trend. Through- 
out this study r=10 was selected as the threshold between tower and 
anvils. 
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Fig. 3. Relationships between deep convective cloud characteristics and aerosol loading. Top row refers to the Atlantic study region (0° N to 
14° N; 18° W to 45° W). Bottom row refers to the Pacific study region (0° N to 14° N; 180° W to 150° W). Left column - (a) and (d), tower 
to anvil area ratio as a function of aerosol optical depth (AOD). Center column - (b) and (e), cloud fraction as a function of AOD. The blue 
curve shows the total cloud fraction. The red shows only the anvil cloud fraction. Note the similarity in slopes. The change of total cloud 
fr action with aerosol is driven by the expansion of the anvils, and the decrease in the relative contribution of the towers. Right column, (c) 
and (f), anvil height and anvil optical depth for the four levels of aerosol loading marked in left figures. The low level AOD is associated 
with die lower right dot. The black arrow marks the AOD increase direction. Increase in aerosol loading is associated with increase of anvil 
height and with a decrease in anvil optical depth or a thinning of the anvil as it spreads out. 


stretches the anvils to cover a larger area, the ice path is di- 
luted in the process, resulting in a smaller TAR (Fig. 3a and 
d) and a decrease in anvil optical depth (Fig. 3c and f). 

2.2 Top of the atmosphere forcing estimations - 
introducing the r-Z cloud forcing space 

How does such a chain of events affect the climate radiative 
forcing? Clouds cool the atmosphere by reflecting back to 
space part of the incoming shortwave solar radiation. They 
warm the atmosphere by absorbing longwave radiation emit- 
ted from the surface and lower atmosphere, and therefore, 
reduce the thermal energy loss to space. These two different 
radiative processes depend on different cloud properties. The 
longwave cloud radiative effect depends on cloud tempera- 
ture and emissivity. Cloud temperature is linked to the height 
and physical thickness of the cloud, and the local profile of 
air temperature. Cloud emissivity is linked to cloud liquid 
water content, but saturates quickly so that cloud emissivity 
is mostly constant for a wide range of cloud optical depths. 
The important point is that higher clouds with colder cloud 
top temperatures will emit less longwave radiation to space. 


thereby warming the atmosphere. In the solar range cloud 
height plays a minor role and the liquid water/ice content, 
cloud optical depth, droplet size distribution and thermody- 
namic phase are the critical cloud parameters that determine 
cloud reflectance of solar radiation. Liquid water content is 
not measured directly from passive remote sensing but can 
be estimated as the product of the cloud optical thickness r 
and cloud effective radius r e , (which is the ratio of the 3rd 
and 2nd moments of the droplet size distribution). Both r 
and r e are retrieved simultaneously from spectral reflectance 
measurements using one visible channel and one channel in 
the mid-infrared (Nakajima and King, 1990: Platnick et al., 
2003). Optically thicker clouds with smaller droplets/ice par- 
ticles will reflect more shortwave radiation back to space. 
The net radiation effect is a superposition of the longwave 
and shortwave processes. 

The different sensitivity of the solar and thermal radia- 
tive regimes to cloud properties implies different response 
in terms of radiative forcing. To demonstrate the sensitivity 
of the top-of-the-atmosphere (TOA) forcing to the deep con- 
vective cloud properties we introduce the r-Z cloud forcing 
space in which each point on this space represents the TOA 


Atmos. Chem. Phys., 10, 5001-5010, 2010 


www. atmos-chem-phys .net/ 10/5001 /20 1 0/ 


I. Koren et al.: Aerosol-induced changes of convective cloud anvils 


5005 


forcing per cloud top height (vertical axis Z) and cloud opti- 
cal depth (horizontal axis r). 

Cloud shortwave reflectance depends not only on the cloud 
optical depth (as presented in the r-Z space) but also on 
the clouds’ droplet effective radius and phase (water or ice). 
Cloud r e will be smaller in polluted conditions due to higher 
concentrations of cloud condensation nuclei (CCN), but will 
increase with height in the cloud. Cloud drop freezing (ther- 
modynamic phase) depends on the temperature, the ice nu- 
clei concentration and the size of the drops. It was shown that 
higher CCN concentration produce smaller water drops and 
by that delay the freezing level to higher altitude and colder 
temperature (Rosenfeld and Woodley, 2000). 

Figure 4a shows the r-Z space of a deep convective trop- 
ical cloud. It shows the daily average net radiative forcing 
calculated using the Santa Barbara DISTORT Atmospheric 
Radiative Transfer Model (SBDART, Ricchiazzi et al., 1998) 
for a typical tropical temperature profile, as a function of the 
cloud optical depth and height. 

It can be seen that any transition towards the left and/or 
top of the r-Z space will be translated as a positive forc- 
ing (wanning). Moreover the gradients in the plane are far 
from being linear, suggesting that the estimation of the aver- 
age forcing from different cloud elements may not be sim- 
ple. To illustrate the above point in a very simple exam- 
ple, Figure 4b shows one cross-section of the r-Z space for 
Z = 10.5 km (marked as a dashed line in Fig. 4a). Note that 
the r scale is logarithmic suggesting that the apparent lin- 
ear slope on the right part of the figure is actually an expo- 
nent. Therefore, since such function is a convex function, 
the Jensen’s inequality applies (Newman et al., 1995). If the 
average cloud optical depth r c in a pixel equals the normal- 
ized by area average of all the clouds in the pixel, the overall 
cloud forcing of this pixel will be larger than the forcmg of 
the average cloud. As an example, say that within a given 
pixel all of the clouds tops are at Z = 10.5 km. Say that a 
cloud with one area element has optical depth z t = 44 and 
another cloud occupies 9 area elements with r a = 4. There- 
fore t c = 8 and from Fig. 4b one can see that the correspond- 
ing forcing in the cloudy area with Z = 10.5 km and r c = 8 
would be F =~ 5 W/m 2 (the magenta line in Fig. 4b). In fact 
the true forcing will be the normalized by area average of the 
forcing namely (9-35-1 -68)/10=24 W/m 2 (see the red lines in 
Fig. 4b) almost 5 times larger than the forcing of the average. 

The sensitivity of the r-Z space on r e is demonstrated in 
Fig. 5. Roughly it shows that the zero forcing line mark- 
ing the transition from negative forcing (cooling, cold colors) 
to positive (wanning, hot colors) occurs higher in the atmo- 
sphere for larger r e , suggesting that for a given r the shift to 
warming will occur m lower clouds in polluted conditions. 
The zero forcing line occurs at lower levels for ice as com- 
pared to water (for a given r and r e ) suggesting that delay in 
freezing will delay the transition to warming. 

For the r-Z space used in Fig. 4, we chose a combination 
of parameters that will mark the higher (conservative) bound- 




Fig. 4. Upper r-Z space of the cloud radiative forcmg at the top 
of the atmosphere, where r denotes the cloud optical depth and Z, 
the cloud top altitude. The color indicate the top of the atmosphere 
(TOA) daily average cloud radiation effect in [W/m 2 ] as a function 
of the cloud height (vertical axis) and cloud optical depth (horizon- 
tal axis) for average effective radius ranging from 20 to 40 pm in the 
upper atmosphere and mixed phase starting from —20 ° to —40 °C 
where full freezing occurs. Warm colors represent a net wanning 
effect, white represents zero net effect and cold colors represent net 
cooling. Any shift of the clouds to the left and/or up in the r-Z 
space is associated with a shift toward wanning. Note how com- 
plex the space is, suggesting that averaging of cloud properties (r 
and Z) will result in wrong radiative flux estimations. Lower: cross 
section taken from the above r-Z space at Z=10.5 km, showing top- 
of-atmosphere (TOA) forcing as a function of cloud optical depth. 
The bottom figure demonstrates that the mean forcing of anvils with 
r=4 (35 Wm -2 ), covering 90% of the scene, and towers with r=44 
(—68 Win -2 ), covering 10% of the scene is ~25 Wm -2 , is 5 times 
higher than the forcing calculated for the mean r of the scene r=8, 
(5 Win -2 ). 


aries from cooling to warming, and will accentuate the cloud 
cooling regime. The transition to ice is delayed until close to 
—40 °C for complete freezing and the effective radius profile 
is skewed towards large droplets and crystals, from 20 pm in 
the lower cloud to 40 pm at higher levels. 
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Fig. 5. Sensitivity of the cloud r-Z space to the cloud thermodynamic phase and the effective radius of the cloud particles (r e ). (a) Water 
clouds with r e = 40 ^im. (b) Water clouds with r e = 20 pm. (c) Only ice particles with r e = 40 pm. (d) Corresponds to Fig. 3 with mixed phase 
and a range of particle sizes from r e = 20 pm at low altitudes to 40 pm at upper altitudes. The white color represents the zero net forcing and 
marks the transition from negative forcing (cooling, cold colors) to positive (wanning, hot colors). Following the location of the transition 
zone reveals that the tendency is for the wanning effect to extend to lower altitudes for ice phase and smaller particles. 


Traditionally, aerosol indirect effects are thought to be pri- 
marily cooling. Here we show that by enhancing anvil devel- 
opment, the result can be significant warming instead. 


3 Discussion 

3.1 Retrieval artifacts and meteorology 

When studying cloud aerosol interactions from observations 
two inherent problems will always surface, namely 1) cloud 
contamination effects on the aerosol retrievals and 2) cause 
and effect - separation of true net aerosol effects on clouds 
from meteorological effects driving both aerosol and clouds. 

On one hand one wants to measure aerosol as close as pos- 
sible to clouds in order to truly reflect the relevant aerosol 
(that interacts with the clouds) properties. On the other hand, 
one asks for very accurate measurements of aerosol load- 
ing and properties, which are very difficult to achieve in the 
vicinity of clouds (Koren et al., 2007, 2008b; Charlson et al., 
2007; Wen et al., 2007). Even if assuming that one can mea- 
sure correctly clouds and aerosol in a cloud field, the next 
challenge is due to the strong coupling of both to meteorol- 
ogy (environmental properties). Are the trends seen between 
aerosol loading and cloud anvils due to real aerosol effects 


on clouds, or does the meteorology drive the changes in both 
aerosol and cloud properties? 

The question of decoupling the meteorology from the net 
aerosol effect was tackled in Koren et al. (2005). There, they 
suggested that vertical velocity from a reanalysis model can 
serve as a good proxy for deep convection. By dividing the 
data to subsets with similar vertical velocity (magnitude and 
sign) and analyzing each set separately they showed that the 
aerosol effect of invigoration of the convective clouds is true 
for any of the subsets. 

A recent study was conducted with the objective of pro- 
viding deeper and detailed answers to these two key ques- 
tions over the same study area (Koren et al., 2010). It shows 
that even after rejection of most of the pixels that may con- 
tain cloud contamination, the trends between aerosol and 
clouds are kept. Moreover, cloud properties as a function 
of aerosol loading show similar trends when using the God- 
dard Chemistry Aerosol Radiation and Transport (GOCART) 
model (Chin et al., 2000) as a measure for the aerosol load- 
ing (instead of satellite measurements). Moving on from 
the question of cloud contamination in the aerosol data set 
to the question of separating meteorology from aerosol ef- 
fects, the study analyzed 280 meteorological variables to see 
which variables best represent the properties of deep con- 
vective clouds measured from space. The selected variables 
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Fig. 6. Same as Fig. 3 except for using two alternative methods to determine aerosol loading. The first and third rows use the MODIS AOD 
data but extrapolate from retrieved grid boxes into fully overcast grid boxes that did not contribute to the statistics of Fig. 3. The first and 
second rows show results of the Atlantic set and the third and forth show results of the Pacific. This method allows inclusion of the tallest, 
most invigorated convective systems. The second and fourth rows use the GOCART chemical transport model (Chin et al., 2000) to provide 
the aerosol information that is collocated with the satellite-retrieved cloud variables. The trends exhibited by these two alternative methods 
are the same as shown in Fig. 3, demonstrating the robustness of the phenomena described. 


helped to divide the data into subsets that represent different 
meteorological conditions. Convective clouds and aerosol 
analysis were conducted per each meteorological state show- 
ing similar aerosol effects for all the sets. 

3.2 Fully overcast pixels 

An additional limitation of this study is related to the ne- 
cessity of having both clouds and aerosol information in the 
same region. As cloud cover approaches full overcast condi- 
tions, aerosol retrievals become rare. Such is the case with 
the most vigorous convection having the broadest spread of 
anvils. This may suggest that our analysis is biased to the 
boundaries of the developed clouds and that we may miss 
part of the more developed clouds that are away from the 
borders. Such limitation may result in underestimation of 
the clouds invigorated by aerosols into the higher atmosphere 
warming regime of r-Z space (Fig. 4a). In order to check 
tills we extrapolated the AOD data to pixels with no data 
as long as they are adjacent to at least one pixel with data. 
Therefore we are estimating the AOD of pixels obscured by 
clouds by averaging the AOD of their neighbour pixels with 
AOD retrievals. Moreover similarly to Koren et al., 2010, the 
anvil analyses was re-done using GOCART aerosol transport 
model information (Chin et al., 2000) as a measure of the 
aerosol loading (instead of the satellite measurements). 


Figure 6 shows both extrapolated and GOCART analysis 
for the Atlantic and the Pacific in a similar manner to Fig. 3. 
The overall trend of the tower to anvil ratio, the cloud frac- 
tions and the cloud height is similar for the extrapolated data 
and nosier for the GOCART. Note that for the Atlantic data 
the TAR is larger for the extrapolated data suggesting that 
the non-extrapolated data underestimate the cloud thickness. 
However, this is not the case for the Pacific data. 

Histograms of cloud fraction (right) and cloud top height 
(left) over the Atlantic study area are plotted in Fig. 7. It 
shows the distributions of all cloudy pixels (red) vs. the ones 
that have both aerosol and cloud information (blue). When 
using the AOD pixels without extrapolation (upper row) the 
analysis misses a significant part of the high clouds with the 
larger cloud fraction. Based on the cloud top height (left) 
histogram it can be seen that most of the effect of the high 
clouds (wanning effect) is likely not to be included in such 
analysis. Wfren using the extrapolated data (lower row) a 
greater portion of pixels with large cloud fraction is included 
increasing significantly the contribution of high clouds. 

3.3 Restricting meteorological variance 

It is suggested in Koren et al. (2010) that the most im- 
portant meteorological variables with respect to spanning 
the variance of deep convective clouds over the Atlantic 
are the broad-scale vertical velocity at mid troposphere (co 
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Fig. 7. Histograms of cloud fraction (right) and cloud top height (left) over the Atlantic study area. Red shows all data and blue shows pixels 
that have both cloud and aerosol information. The upper row shows data used for the analysis using the MODIS AOD without extrapolation 
(Fig. 3) and the lower row shows data with extrapolation (Fig. 6). It shows how the extrapolation (done only in cases where the adjacent 
neighbour has a valid AOD retrieval) of overcast pixels adds significantly to increasing the number of high clouds in the analysis (Fig. 7c). 



Fig. 8. Same analysis as shown in Fig. 3, but dividing the data into two meteorological regimes defined as either large scale rising motion 
(negative pressure vertical velocity, co ) and denoted as solid cmves, or subsidence (positive co) shown as dotted curves. 


550 hPa) and the relative humidity in the upper atmosphere 
(RH 350 hPa). Based on this, here we restrict meteorological 
variance using the mid-tropospheric vertical velocity from 
the NOAA-NCEP Global Data Assimilation System (GDAS, 
Kanamitsu, 1989: Parrish and Derber, 1992). Figure 8 shows 
the analysis of Atlantic data similarly to Fig. 3, but the data 
is divided into two sets with positive (dotted) and negative 
vertical velocity. Note that co is defined as a pressure verti- 
cal velocity so that negative values imply rising motion and 
positive values denote subsidence. Figure 8 shows that in- 
deed the clouds with negative co are more developed, taller 
and cover larger areas. It is also shown that the aerosol effect 
is ahnost orthogonal to the meteorological subsetting. The 
same trends exist between the TAR and cloud fraction with 
AOD, and that the effect of the meteorology is a simple offset 
of the trend line. The trend in anvil optical depth (Fig. 8c) is 
less profound for the positive co set. 


4 Summary 

In numerous climate studies cloud aerosol interaction is 
shown to have the potential to be a key contributor to the 
radiative forcing (see Heintzenberg and Charlson (2009) for 
detailed information). Small changes in the cloud shape, 
structure or lifetime can change significantly the local radia- 
tive balance. 

Here we suggest that since the wind speed increases sig- 
nificantly in the upper atmosphere, if indeed aerosols invig- 
orate deep convective clouds, such invigoration will have the 
following consequences: invigorated convection associated 
with increased aerosol loading results in expanded anvil spa- 
tial coverage and reduced anvil optical depth in higher levels 
of the atmosphere. These thinner, colder clouds might induce 
strong climate wanning. 
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We introduced the r-Z space of cloud forcing and showed 
that such forcing representation reveals the opposing effects 
and the different sensitivity on cloud properties of the so- 
lar vs. the thermal spectral regimes. Moving left in the di- 
agram - towards lower r, or up - towards higher altitudes 
indicates transition to warming. The aerosol effect suggested 
in this paper is moving clouds both left and up in the r-Z 
space. Moreover, we have shown a more general problem 
in cloud forcing estimations namely since the forcing depen- 
dence on both the cloud optical depth (r) and height (Z) is a 
convex function in a significant portion of the space, the forc- 
ing of the averages (r or Z) is much smaller than the average 
forcing suggesting forcing underestimation at any scale other 
than that of the individual cloud type. In particular the r-Z 
space suggests that the anvil and the tower have to be treated 
independently for the forcing calculations. 

The limitations of such analysis were discussed, hi short, 
we try to study aerosol interaction on deep convective clouds 
while knowing that measuring aerosol properties in the vicin- 
ity of such clouds is challenging. We devote a whole study 
specifically to examine such problems (Koren et al., 2010). 
Here we find the same associations between aerosol and con- 
vective cloud characteristics in two different geographical lo- 
cations with entirely different aerosol transport and different 
oceanic surface characteristics. We find these same associ- 
ations using two different methods to collocate aerosol and 
cloud properties and when we restrict meteorological vari- 
ance. We even find the same relationships between aerosol 
and convective cloud properties when we use a model to 
provide the aerosol loading information. The results of this 
study are reproducible and robust, and suggest strongly that 
the proposed trend is a true aerosol effect. 

Due to data limitations however, we do not quantify the 
radiative forcing of the proposed mechanism. Though crit- 
ically important, reliable estimation of the radiative forcing 
due to aerosol effects on deep convective clouds and their 
anvils will have to be done by using a more complete set of 
information on clouds, aerosol and meteorology, in conjunc- 
tion with models. Here we demonstrate a pathway of aerosol 
effects on clouds that can lead to a positive climate forcing. 
This is an important mechanism that will be required in fu- 
ture estimates of climate change. 
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